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Abstract 

We analyze a class of stochastically stable quenched measures. We prove that stochas- 
tic stability is fully characterized by an infinite family of zero average polynomials in the 
covariance matrix entries. 
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1. Introduction 



After three decades from their first appearance in the Edwards and Anderson work [EAJ 
the spin glass models and their low temperature phase remain one of the major unsolved 
problems of condensed matter physics. The physically well understood mean field case (see 
Parisi et al. [MPV]) is still under investigation from the mathematically rigorous perspective 
and some recent results by Guerra and by Talagrand [Gl' TS] confirm the Parisi theory. The 
case of spin glasses on finite dimensional lattices is instead much more controversial, and the 
structure of its equilibrium states is an unsettled matter even from the theoretical physics 
point of view. 

The spin glass model problems arise from a peculiar mathematical structure of two 
intertwined probability measures, the configurational (spins) and the disordered (random 
couplings) which are combined in a precise measure prescription of equilibrium statistical 
mechanics commonly called quenched ensemble. 

In Aizenman and Contucci [AC] a stability property for the mean field models was 
derived from the continuity (in the temperature) of the thermodynamic functions. Some of 
the consequences of such a stability were proved to be captured by an infinite family of zero 
average polynomials (see also Ghirlanda and Guerra [GGJ). Subsequently in jT] the same 
property was investigated in finite dimensional models and proved to imply a formally similar 
property in terms not of the standard overlap function but of the so called link-overlap. 

The stability property, nowadays called stochastic stability, attracted some attention 
from both theoretical and mathematical physics. It was first investigated in Franz et al. 
|FMPP1I [FVIPP 2J and cleverly used to determine a relation between the off-equilibrium 
dynamics and the static properties. More recently a purely probabilistic version of it expressed 
in terms of invariance under reshuffling of random measure for points in the real line has 
been investigated and completely classified in the case of independent jump distribution (see 
Ruzmaikina and Aizenman [RAJ). The connection between the probabilistic approach and 
the statistical mechanics one is well explained in | |G2j and based on a new variational principle 
introduced in [ASSJ. 

Can we give a complete characterization of the stochastic stability property within its 
original statistical mechanics formulation? In other terms once we know that a spin glass 
model verifies stochastic stability do we know what are (all) the constraints of its overlap 
distribution? In this paper we answer positively the previous questions and we prove that 
thanks to a remarkable cancellation mechanics already observed in jC2| . the zero overlap 
polynomials of [AC] or [CJ provide a complete description of the mentioned stability property. 

The paper is organized as follows: in Section 2 the quenched measure is introduced and 
the overlap moments formalism explained. Section 3 introduces a combinatorial description 
of the overlap measure on graph theoretical grounds. Section 4 introduces stochastic stability 
and contains the main result (Theorem [TU]). It states a property which implies that all the 
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consequences of the stochastic stability are indeed contained in its second order version. 



2. Quenched measures 

A quenched probability space is a product measurable space rijx^l^, where the random 
probability measure fij on is indexed by J G ^Ij, and distributed according to a probability 
measure v on Qj. 

Example 1 The Sherrington-Kirkpatrick (SK) model with spins (no external field). 
Here Q„ := { — 1,1}^, Qj := M^^, random interactions Jij with Cartesian coordinates 
1 < ^, j < ^ and Gaussian measure u of density (27r)~^ ^'^ J2i<i j<N '^Ij /'^) w.r.t. 

Lebesgue measure. For inverse temperature P > the random Gibbs measure J i-^ /xj is 
given by 

:= (- - • • • ' e n.), (2.1) 

Hjia):=- J2 AM^Hj)^ (2-2) 

^<i,j<N' 

and partition function Zj{/3) := J2aen^^^P(~(^'^j('^))- 

Example 2 The Edwards-Anderson (EA) model with spins in a volume A C Z*^ and 
nearest neighbor interactions. 

Indicating by 5(A) the set of nearest neighbors of A: := ilj := RI-^^^)!, 

random interactions J(j j) for nearest neighbors in Cartesian coordinates E B{A) and 
Gaussian measure u of density (27r)"l^('^)l/^ exp(— ^.-j^^^^^ J!^^/2) w.r.t. Lebesgue mea- 
sure. As in the former example the random Gibbs measure is given by ()2.ip with the Hamil- 
tonian 

Hjia):=- J2 J-^M^^U)- (2.3) 

(ij)GB(A) 

For random variables / : i7j x i^^. ^ R resp. (7 : llj — > R we use the notation 



{f){J) ■■= / /(a,J)ci/ij(a) , Avig) := / giJ) di^iJ) . 

Quantities of particular interest in a quenched probability space are the moments ({fY) 
(r G N). We denote the elements of the product space of i? G N real replica by 
cr = ((Ti, . . . , aji), and equip it with product random measure J ^ ji^ = ® . . . ® fij. 
Expectation w.r.t. fij is denoted hy <^ — ^ (J). 
The (non-random) quenched measure on ilj x has expectation 

E(-) := Av(< - >). 
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It is not necessary to specify the number R of replica, as a function on Vl^ can be extended 
(by inverse projection) to 

Our plan is to study a specific class of quenched models whose algebra of observables 
is built over families of Gaussian variables. Later on we additionally assume that the two 
measures involved fulfill a remarkable stability property. 

By enlarging (f2j,z/) we introduce additional standard normal random variables 

h^''\a) -.Qj indexed by (7 G and A; = 1,. . (2.4) 

We assume that the h^^\a) are independent of the former J variables, that h^^\a) is 
independent from h^^'\a') \^ k ^ k' and and that their covariances 

do not depend on A; G {1, ... , K}. 



Example 3 For the SK model of Example[T]and independent standard normal variables J, 
(i, j = 1, . . . , A^; /c = 1, . . . , K) the Gaussian variables 



(fc) 



N 



have covariances Av(^h^''\a)h^'''\a')^ = 5k k'Caa' with 2^ x 2^ matrix entries explicitly given 

by 

c.,.' := (^-^Yl^a{i)a'{i)^ G [0,1]. 
Example 4 Equivalently for the EA model of Example Elthe Gaussian variables 

have covariances 

c^y := \B{A)r Yl <i)<jy{i)^\3) ^ ["1, !]• 
(i,i)eB(A) 

Finally, for the replica space Vt^ we introduce 

Cr,r' ■ ^ [^l^ 1] ' Cr,r' (c^l, . . . , (Tr) := C^^^, (r, r' = 1, . . . , i?) . 
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The indices r enumerate the replica. By the normality assumption 

Cr,r{cr) = 1 {r = 1, . . . ,R; a e Q^). 



Using Wick's Theorem for a family of Gaussian random variables, see Glimm and Jaffe 
|GJ| or Simon 0, expectations of products of h^^\a) lead to sums of products of Cry: 
Whereas averages over odd products vanish, for R = 2m 

(R \ m 

1=1 / pairings tt i=l 

m m 

pairings tt i=l pairings vr i=l 

Here a permutation tx -.V -^V \s called a pairing \N\th ordered pairs 

(7r(l),7r(2)),...,(7r(i?-l),7r(i?)) 
if 7i{2i - 1) < 7r(2i) and n{2i - 1) < 7r(2z + 1). There are {R - 1)!! pairings of V. 
Example 5 



Av{{hf)= [ du{J)\( [ diij{a^)h{a,,J)\( f dfij{a2)h{a2, J)] 



E(ci_ 



2J, 



and 

Av((/i«>(/i«/i(2))(/,(2))) = 

rfz/(J)f / rf/ij(ai)c?/ij(a2)rf/xj(a3)/i(^)(ai, J)h^^\a2, J)h^'\a2, J)h^'\as, J) 

= E(ci,2 C2,3) • 

The general task will thus be to analyze the expectations of the so-called overlap 
monomials, that is, random variables on the replica space fi^ of the form 

Yl c^^'' with nii^j = nij^, G Nq. (2.5) 

l<j<j<-R 

For any permutation tt : \^ ^ ^ of the replica index set V := {1, . . . , R} 

e( II cZ^^^^-A = e( H (2.6) 

\i<i<j<R / Vl<i<j<K / 



(e.g. lE(ci,2C2, 301,4) = IE(ci,2C2,3C3,4)), since the random probability measure fXj is vr-invariant. 

The main observation of this section is that the property of the quenched probability 
space can be studied as the probability expectation E of the random covariance matrix C. 

Next section introduces the suitable combinatorial language for the study of the overlap 
distribution in terms of its moments. 



3. Multigraphs and the Gaussian Operator 5 

In order to analyze expectations of the overlap monomials ()2.5p . it is advisable to use 
the notion of (edge-labeled) multigraphs (see e.g. Diestel [D] for more information). 

Definition 6 A multigraph (on R eN vertices) is a finite set E and a map 

G:E^ [V]^ U V 

with the vertex set V = {1, . . . , R} and the family [V]"^ := {{i,j} \ i 7^ j ^ V} of unordered 
vertex pairs. 

We call the elements of E' := G-\V) legs, the elements of E" := G-\[V]'^) = E - E' 
edges. 

The set of all multigraphs is denoted by Q, and 

00 

g= [j with := {G I \E"\ = m, \E'\ = n}. 

m,n=0 

We write the multiplicities [(^^^(e)! as exponents: 

Example 7 G = {1, 2}2{1, 3}{2} G ^^'^ is the graph with the edge {1,2} of multiplicity 
2, the edge {1,3} and the leg {2} of multiplicity 1, and G G G^'^'^l 

We use the shorthand Q" := [J.^G^"^'^^ for the subset of multigraphs without legs, and 
similarly Q' := IJn multigraphs without edges. 

Q is the basis of the vector space Q := C[Q] of finite linear combinations, and we 
denote subspaces spanned by the various subsets of ^ by a tilde (using multiplication of 
multigraphs and formal linear combinations, Q becomes an algebra, sometimes called overlap 
algebra) . 

We now introduce the linear operator C : Q ^ Q of Wick contraction (the name 
being justified by (13. 2 1) below) by 

m 

C{E) := E"C{E') , C{E') := H^^^^^^^ - 1)), G(7r(2z))} . (3.1) 

pairings tt i=l 
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for the decomposition E = E' U E" of the multigraph G E Q (here for \E'\ = 2m the 
pairings are seen as permutations tt : {1, . . . , 2m} — > £" = {1, . . . , 2m}). Note that 

C = {0} and C c ^(-+'.0) . 

Using the bijection I : G ^ G between overlap monomials 

n ^•]{h>A 

\i<i<j<R / \i=l / 

and multigraphs G : E [V]'^ U V with edge multiplicities \G^^{i,j)\ = mij and leg 
multiplicities \G^^{i)\ = rii , vje see that 

IAv{G) =CI{G). (3.2) 

Thus we omit the hat of G in the rest or the article. The invariance (12. 6 p under permutations 
allows us to freely use the multigraph isomorphisms induced by arbitrary relabeling of the 
vertices in order to calculate expectations. 

In order to prepare for the deformation of the random measure treated in the next 
section, we introduce a second linear operator: 

6 : g^g , 5G := J2 ^-^ ^ 

veViG) 

with 

6, := 5W + , and ;= {v}G , di'^G := -{v'}G, (3.3) 

where v' is the first element of N not belonging to F (so t>' = i? + 1 if ^ = {1, . . . , R}). 
Note that 

and the properties 

5iGiG2) = {5Gi)G2 + Gi5G2 , 5(0) = 0. (3.4) 

of a derivation (0 denoting the graph without legs and edges, the neutral element under 
composition of multigraphs). 

Example 8 • <5^^^{1, 2} = {2}{1, 2} 
. 4-){l,3} = -{2}{l,3} 
. 4-)4-){l,3} = {4}{2}{l,3}. 
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4. Stable Quenched Measures 



In this section we investigate the consequences of a stability assumption for the couple 
(E, C) for suitable random deformation. Stochastic stability holds, strictly speaking, for mean 
field and finite dimensional models once the thermodynamic limit is taken. While studying 
its consequences we will have to consider derivatives of the moments of (E, C) with respect 
to suitable parameters: our computations will tacitly use the fact that the thermodynamic 
limit does commute with the operation of derivative as it was proved, for every order of 
derivation, in the appendix of [ACJ . 

Given a quenched probability space with a family ()2.4p of standard Gaussian random 
variables we introduce the deformed random measure with expectations 

and the corresponding deformed quenched measure: 

ExhH ■■= Av(< - :^xh) , (4.2) 

where h = h{a, J) is a Gaussian random variable of covariance c^y, //-independent from 
the family (h^-^\ . . . ,h^^^), say h = h^^^^\ and A G M parameterizes the deviation from 

E = Eq. 

We denote a monomial in the covariance matrix entries by G. Its expectation is always 
an even function of the parameter: 

E_,,(G)=E,,(G') (AgM). (4.3) 



The quenched probability space defined in the former section is defined to be stable 
when (E, C) is isomorphic to {¥.\h,C) for every A or, since the entries of C are bounded, 
when the moments of the two measures E and E^h coincide for every A. Stability of the 
measure E under deformation means that for every monomial G 

Exh{G) = E(G') , (4.4) 

and in particular, as the left hand side is A-independent: for all G 



Q2n 

-E;,,(G)U=o = 



(we consider only the even derivatives because the odd ones vanish independently of the 
assumption (14. 4p by the symmetry of the Gaussian). 

Whereas typically finite spin systems are not stable, for systems like the SK model 
stability is thought of being an asymptotic property in the thermodynamic limit — > oo. 
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Proposition 9 For every G E Q" and every integer n 

rEA,(G)U=o = E(C52"G) 



QX^n 

where C5^"G is a polynomial in the matrix entries. So the stability property (14. 4|) implies 

E(C5^"G) = 
for every G and every n and in particular, defining 

A := C5^ 

that _ 

E(AG') = {Geg"). (4.5) 

Proof: (See also Prop. 5.1 in |AC| ) The action of the operator 5 on the multigraph basis 
corresponds to the usual derivative with respect to the parameter A: 



d 

dX 

So such a derivative produces a truncated correlation expressed in the rule ()3.3p . 
The derivation property ()3.4|) represents the Leibniz rule 

R 



d 
dX 



«-»- = e(^hS)(8)(-) 

1=1 ^ ^ r^l 



(r) 
\h 



for the derivative on the replica space. Each differentiation with respect to A produces a sum 

of monomials, each containing an additional zero mean Gaussian variable h. 

Finally, the contraction C represents Wick's Theorem. □ 

Our main result is 
Theorem 10 

= (2n-l)!!A" (neNo), (4.6) 

or equivalently 

Q2n _ 



-E,,(G)U=o = (2n - 1)!!E(A"G) {n G No, G G G") . (4.7) 



As a consequence the stability property is equivalent to the set of relations (14. 5|) 
Remark 11 In the sense of formal power series we can combine the set of identities (l4 



as 



Cexp(t6) = exp(iAt2) (t G C), (4.8) 
which justifies to call S Gaussian operator. 
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The proof is based on the following: 
Lemma 12 For all G e Q" , we find (with A = C5'^) 

-^^xhiG) = mni^G) (4.9) 
Proof: By definition of 5, we know that for the overlap monomial the l.h.s. equals 

^^xh{G) = Eif,{6G) =Eih (^{-Rhjt+, + J2h,)G^ , (4.10) 
and that the r.h.s. of ()4.9|) equals 

Exh{AG) = Eih ( ( 2 5^ Cij - 2i?^ + R{R + 1))ch+i,r+2 ) G J , (4.11) 

V V l<i<j<R i=l ) ) 

see |AC] . Lemma 6.2. 

We now do integration by parts in order to get rid of the unpaired Gaussian variable 
that appears in the r.h.s. of (I4.10|) . The same integration by parts is responsible for the 
appearence of the factor A in ()4.9|) . Recall that for a family (h^^\ . . . , h^^^) of Gaussian 
variables with covariance matrix {ci^m} the following rule of integration by parts 

e(/.«/(/.«,...,/.(^))) = X:e(q,„^/(/.W,...,/.W)) (4.12) 

m=l ^ 

holds (if / G C^(]R'^) is of moderate growth at infinity like for h i-^ exp(A/i)/ (exp(A/i))), 
see e.g. [CDGGJ or Talagrand [TlJ, A. 6. 

Applying PTT^ to I^TW^ we obtain I^TT^i . □ 



Example 13 Let us explicitly verify the lemma in the easiest non-trivial case G = {1,2} 
(see also Sec. 5 in [ACJ and Theorem 2 in [GJ). From (14.111) we see that 

C6'{1, 2} = 2 [{1, 2}2 - 4{1, 2}{2, 3} + 3{1, 2}{3, 4}] . 
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Applying to G the first derivative in A and then integration by parts formula ()4.12|) . we find: 



d_ 

dX 



^xhicu) = Av 



Cuih + h2) 2 / Ci2/l3- 



2A Av 



Cl2C2j 



Q /gA(/li+ft.2) 



Cl2 ( (C21 + 1) 



gA(/li+/l2) 



Cl2C3j 



Z2 



gA(/li+/l2+/l3) 

2C23- 



Z3 

^ / ^\{hi+h2+h-i) 

dhi 



gA{/li+/l2+/l3) 

Cl2 ( (Csi + C32 + 1) 3C34 



'^l.°"2.°"3''^4 



^'X(h-i+h2+hz+hi) 



Z3 



Z4 



2A Ea/i [Ci2 - 4ci2Ci3 + 3C12C34] 



The associated multigraph polynomial of this expression, 2[{1, 2}^— 4{1, 2}{2, 3}+3{l, 2}{3, 4}], 
is just equal to C5'^G. 

Proof of Theorem llOt Considering G E G" and thanks to equality ()4.9p . we can write the 
next derivatives in A as follows: 



and for general /c G N by induction 

^E,,(G) = (fc-l)^E,,(AG) + /^E,.(AG). 
Considering that in / = all the odd derivatives vanish by ()4.3p . we obtain: 



& 



2n 



2n-2 



E,,(G)b=o = (2n-l)— — E,,(AG)|,=o. 



Combining eqs. (I4.15p 



2n 



E,,(G)| 



gi2n-2 



(2n- 1)!!E(A"G). 



(4.13) 



(4.14) 



(4.15) 



In terms of the operation 5 the previous relation becomes 

C5'" = (2n - 1)C5'"-2A = ... = (2n - 1)!! A" , 

which proves (14. 7j) . 

To prove equivalence of ()4.5jl with stability, we observe that A"G is a polynomial in the 
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overlap algebra i.e. A^G = J2a ^G^C^a so that the vanishing of the average of AG for every 
G implies also the vanishing of the average of A"G: 

^E,(G)U=o = (2n-l)!!E(A"G) = (2n-l)!!5^nc„E(G„) = (4.16) 

a 

Remark 14 In order to fully appreciate the content of the ()4.6p we examine the case n = 2 
(see also |C2|): 

CS^ = 3A2. (4.17) 

It is interesting to note that the left hand side of (I4.17|) contains a priori 3 ■ 2^ = 48 terms 
(the factor 3 coming from the Wick contraction of a fourth order monomial and 2^ being 
the number of terms of (5+ + S~y ). The right hand side instead contains only 4^ = 16 
terms coming from the square of A. 

Although by definition the Wick contraction does not conserve the number of edges 
nor the number of vertices, the presence of alternating signs in the definition of 5 together 
with the invariance under permutation of graph labellings produces a delicate cancellation 
mechanism responsible for the clean equality ()4.17p . 

Acknowledgments. One of us (PC) thanks M.Aizenman, S. Graffi and F.Guerra for 
many interesting observations on the property of stochastic stability. 
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